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ABSTRACT 

The correlation between the large-scale distribution of galaxies and their spectroscopic 
properties at z = 1.5 is investigated using the Horizon MareNostrum cosmological run. 

We have extracted a large sample of 10 5 galaxies from this large hydrodynamical 
simulation featuring standard galaxy formation physics. Spectral synthesis is applied 
to these single stellar populations to generate spectra and colours for all galaxies. We 
use the skeleton as a tracer of the cosmic web and study how our galaxy catalogue 
depends on the distance to the skeleton. We show that galaxies closer to the skeleton 
tend to be redder, but that the effect is mostly due to the proximity of large haloes 
at the nodes of the skeleton, rather than the filaments themselves. 

This effects translate into a bimodality in the colour distribution of our sample. 
The origin of this bimodality is investigated and seems to follow from the ram pressure 
stripping of satellite galaxies within the more massive clusters of the simulation. 

The virtual catalogues (spectroscopical properties of the MareNostrum galax- 
ies at various redshifts) are available online at http://www.iap.fr/users/pichon/ 
MareNostrum/ catalogues. 

Key words: large-scale structure of the Universe galaxies: evolution methods: N- 
body simulations hydrodynamics 



1 INTRODUCTION 

During the past decade, the ACDM cosmological model of 
the Universe has been established as the framework of choice 
in which to interpret how and when observed galaxies ac- 
quire their properties. Arguably the most important feature 
of this framework is to provide us with an explanation as 
to why many of these properties (physical sizes, luminosi- 
ties) strongly correlate with galaxy mass while others (star 
formation rates, morphological type) do not seem to. Un- 
surprizingly, the all-time favored culprit is the interplay be- 
tween galaxies and the intergalactic medium (IGM) at large. 
In other words, the large scale environment of galaxies is 
claimed to play an important role in shaping some of their 
properties, while the rest of them are thought to depend 
solely on small scale (internal) processes. However, having 
said that, one still has to determine for which of these prop- 
erties "Nurture" dominates over "Nature" and therein lies 
the whole difficulty of the issue. 

Indeed, since the early 70s, there has been a plethora of 
studies devoted to measuring the impact of environment on 
galaxy properties. |Davis fc Geller| ( |1976[ ) first pointed out 
that early-type galaxies are more strongly clustered than 



the late types, while Dressier ( 1980 1 and Postman & Geller 
||1984[) demonstrated the existence of a morphology-density 



relation (MDR). Following in their footsteps, Balogh et al. 


(19981, Hashimoto et al. 


( 1998 1 and more recently Christlein 


& Zabludoff (20051 and 


Poggianti & von der Linden ( 2006 1 



systematically showed that galaxies living in denser envi- 
ronments tend to be redder and have lower star formation 
rates (SFRs) than their more isolated counterparts. One can 
think of several physical processes associated with differ- 
ent types of environment that could play a role in causing 
such alterations. More specifically, they include, by ascend- 
ing order of environment density: (i) major (wet) mergers, 



which can turn spiral galaxies into ellipticals (e.g. Toomre & 



Toomre ( 1972 l) , and drive a massive starburst wind which 



quenches future star formation by ejecting the interstellar 



medium (ISM) out of galaxies (Mihos & Hernquist (19961; 
Mac Low & Ferrara ( 1999 l); (ii) active galactic nuclei (AGN) 
or shock-driven winds (Murray et al. (20051; Springel et al. 
( 12005} ; (hi) galaxy "harassment" (rapid encounters or fiy- 
bys which dominate over mergers in rich clusters) causing 
discs to heat and possibly triggerring the build-up of a bulge 



via the formation of a bar (e.g. Moore et al. ( 1998 1). For the 
latter category, the diffuse gas associated with the galaxies' 



2 C. Gay, C. Pichon, D. Le Borgne, R. Teyssier, T. Sousbie, J. Devriendt 




Figure 1. A 3D view of the galaxies and the skeleton of the dark 
matter of MareNostrum at Z = 1.6. The box is 50 h _1 Mpc aside. 



host dark matter (sub)halo which constitutes the main fuel 
supply for future star formation can also be stripped, thus 
suppressing later star formation by "starvation" or "stran- 
gulation" ( |Larson et al.|1980| |Bekki et al.||2002[ ). Moreover, 
part of their ISM can also be pulled out of these galaxies, 
either by tidal forces arising from the gravitational potential 
of the cluster or by ram pressure stripping by the intraclus- 



ter medium ( I CM)|Gunn fc Gott|fl972| >; |Abadi et al.| ( |1999[ ); 
|Chung et alT] ( |2007[ ). 

Although these latter environment dependent processes 
seem potent enough, recent work carried out by |Tanaka et al.| 
( [2004] ) and van den Bosc h et al.| ([2008 ) indicate that they 
might not be the main mechanisms for quenching star forma- 
tion activity. This claim is corroborated by the higher red- 
shift results (z~ 1) obtained with the DEEP2 (e.g. |Cooper 



et al 



(2008 



iTasca 



(2006 2007 20081; Gerke et al. (20071; Coil et al 



I), (z- ) COSMOS ( Scoville||2007| |Cassata et al.||2007 
2009] ) and VVDS |Scodeggio fe Vergani|20"09f lurveys 



where clusters are more scarce, along with the fact that mor- 
phological and spectrophotometric properties of local galax- 
ies are also found to be correlated with their internal proper- 



ties, such as luminosity, mass or internal velocity (e.g. Kauff- 
mann et al. (20031). Here as well, one can invoke various 



physical processes to explain such dependences on internal 
properties. Supernova feedback can heat and stir the ISM, 
possibly ejecting large amounts of gas out of galaxies and it 
is expected to scale with galaxy mass ( |Larson|1974||Dekel fc| 
Silk|1986[ ). There also exists a growing host of observational 
evidence that AGN play a key part in quenching star forma- 
tion ( |Schawinski1[2056{ |Schawinski et al.||2007| |Salim|[2007l ) 
and this AGN feedback should likely impact more massive 



galaxies since they host larger mass black holes (Magorrian 
|et al.|1998||Silk fc Rees|1998[ ). In light of these investigations, 
it becomes apparent that disentangling nature and nurture 



is a more complicated process than one would naively have 
thought to begin with. 

Clearly, the fundamental requirement to tackle this is- 
sue is to properly characterize the anisotropic environment 
of galaxies, both in observational samples and theoretical 
models, spanning as broad a range of environments as pos- 
sible, from isolated field galaxies to groups and rich clusters. 
The vast majority of the studies in the literature accomplish 
this task either by counting the number of neighbours that 
a galaxy has within a fixed aperture on the sky or by mea- 
suring the distance to the n th nearest galaxy, where n is 
an integer in the range 3-10. Although these indicators are 
straightforward to obtain, their physical interpretation, let 
alone their comparison to theoretical models are far from 



being straightforward (c/. Kauffmann et al. (20041; Wein- 
|mann et a.l.| ( |2006| | ) . Meanwhile, looking at the distribution 
of observed galaxies in modern cosmological surveys, such 
as the 2dF (|Colless||2001|) and the SDSS (|York||2000|), the 



most striking feature is that they look organised along linear 
structures linking clusters together (see Figure [l I . This fila- 



mentary network, dubbed as the "Cosmic Web" ( Bond et al. 
|1996[ ) , has a dynamical origin and reflects the anisotropic ac- 



cretion taking place in clusters (Sousbie et al. 2008 1 . It there- 



fore seems natural to describe the environment of galaxies 
in terms of their location with respect to these filaments in 
order to investigate the influence of the Cosmic Web on the 
properties of the galaxies it encompasses. 

In this paper, we carry out such a study at intermedi- 
ate redshift (z ~ 1.5), mainly from a theoretical perspec- 
tive, using a recent diagnostic tool to characterize the 3D 
environment called the skeleton ( Sousbie et al.|2008[ ), which 
we combine with the largest hydrodynamical cosmological 



simulation performed to date (Ocvirk et al. 2008 Prunet 



|et al.||2008"l |Dekel et aLl|2009| ). Run on the MareNostrum 
computer at the Barcelona Supercomputer Center using the 
RAMSES code (iTeyssierl 120021), this simulation is one of 



the flagship simulations realized by the Horizon collabo- 
ration (http://www.projet-horizon.fr). It includes a de- 
tailed treatment of metal-dependent gas cooling, UV heat- 
ing, star formation, supernovae feedback and metal enrich- 
ment. 

Specifically, we will address the question: are the phys- 
ical conditions within the filaments dramatic enough to 
strongly influence the properties of the galaxies it encom- 
passes? 

The outline of this paper is as follows: first we describe 
in Section [2] our methodology, in terms of numerical tech- 
niques, estimators and statistical measurements. The depen- 
dence of the spectroscopic properties on the filamentary en- 
vironment is then discussed in Section [3] while Section [4] 
investigates the observed bimodality and discusses compar- 
ison to observations, and Section [5] wraps up. Some checks 
are performed in Appendix A, Appendix B describes the 
publicly available catalogues and Appendix C sums up the 
subgrid physics used. 



2 METHODOLOGY 

Let us first describe the MareNostrum simulation, a cosmo- 
logical N body and hydrodynamical simulation of unprece- 
dented scale which accounts for most of the physical pro- 
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Figure 2. a sample of MareNostrum galaxies 
the IRAC-8/xm, K, I bands which were used 
this investigation. More such images are available 
http : //www. iap . f r/users/pichon/galaxies/. 



in 

at 



cesses involved in galaxy formation theory. The procedure 
to identify galaxies and assign spectra and colours to them is 
also described. Finally the method used to identify filaments 
within the simulation is presented. 



2.1 The MareNostrum simulation 

We use the MareNostrum simulation, described in detail in 



Ocvirk et al. (20081. To summarise, this simulation uses the 



AMR code RAMSES |Teyssier| ([2002]) in a periodic box of 
comovmg 50 h' 1 Mpc, with a ACDM universe (Qm = 0.3, 
Q A = 0.7, D.B = 0.045, H = 70 km-s _1 -Mpc _1 , cr 8 = 
0.9). The initial grid consists of 1024 3 dark-matter parti- 
cles (m part . ~ 8 x 10 6 Mq) and the same number of cells, 
which are refined up to five folds when they have more than 
8 particles, as long as the the minimum cell size is not un- 
der 1 kpc in physical units. In addition to its large volume, 
the MareNostrum simulation provides the basic physical in- 
gredients relevant to galaxy formation. The hydrodynamic 
physics uses metal-dependent cooling, UV heating ( jHaardt 
& Madau|(|1996|) background model), star formation ( jRasera 



& Teyssier 200 6}, superno vae feedback and metal enrich- 
ment ( |Dubois fc Teyssier||2008[ ). The ISM, i.e. gas with 
density above 0.1 atom/cm' 5 is modelled with a polytropic 
equation (Schaye & Dalla Vecchia 20071 and forms stars 



consistently with the Kennicutt law with a star formation 
efficiency of 5% (see Appendix [c]). For each snapshot, the 
dark-matter substructures are detected with the Adaptahop 
algorithm ( |Aubert et aT]|20u4| |Tweed et aT1|2009| ) and all 
stars are associated to a (virtual) galaxy. The simulation 
stopped at redshift z w 1.5, with 5 x 10 9 cells, 1.4 x 10 8 
star particles and around 100,000 galaxies. The simulation 
parameters (Lbox = 50 /t -1 Mpc, 1024 3 dark matter parti- 
cles, and a spatial resolution close to 1 /t -1 kpc physical) are 
optimal to capture the most important spectral properties 
of typical Milky- Way-like galaxies. The box size allows us 
to have a large sample of about 100,000 galaxies at red- 
shift 2 and above, with a strong statistical significance (see 
Appendix [B] for a detailed description of our sample). 



2.2 Spectral synthesis 

One main interest of the MareNostrum simulation lies in its 
ability to yield realistic virtual observations in a consistent 
cosmological framework, which in turn can be compared to 
real data. The outcome of these comparisons should lead to 
clues about the physics which drives the evolution of galax- 
ies. 

To make such predictions, light needs to be added to 
the simulation. This can be done very naturally by asso- 
ciating spectral synthesis population models to the star- 
formation modelling in the MareNostrum simulation. More 
precisely, the ISM gas produces a single stellar population 
(SSP) per gas cell at the rate described in section 



2.1 



Each 

of these SSPs has the metallicity of the gas which gave birth 
to the stars. It is described by a particle in the simulation 
with a mass M* (minimum value ~ 1O 6 M0), a metallic- 
ity, a redshift of formation and a position. Then, we assign 
a dust-free evolving spectral energy distribution (SED) to 
each of these stellar particles with the PEGASE.2 ( |Fioc fc 
Rocca-Volmerangc 1997 1999) population synthesis code. A 
Salpeter IMF is used for the spectral modelling, consistently 
with the SN feedback. 

Finally, the light content of a galaxy at redshift z is 
the sum of the SEDs produced by its "stellar" (i.e. SSP) 
particles which depend on their age and metallicity: 

n 

F A (A, *) = ^M*x Fl M ° SSP (A, t(z) - t(4° r ), Z(4° r )) , (1) 

i=0 

where A is the wavelength, n is the number of stellar particles 
inside the galaxy, M* is the total stellar mass of the i th SSP, 
t{z) is the Hubble time at redshift z and zf or is the redshift 
of formation of the i l SSP, and Z is the metallicity of the 
gas. Here, in contrast to the usual approximation of instan- 
taneous mixing of chemical elements in a galaxy adopted 
by many models of galaxy formation (such as semi-analytic 
models, e.g. |Hatton et ak| ( |2003| ), the simulation makes it 
possible to account for spatial variations of the metallicity 
within a galaxy. Recall that in this paper a galaxy is defined 
to be a set of at least 10 stars which are embedded within a 



given dark halo subclump. Rimes (2009) presents an alterna- 



tive definition of a galaxy within the MareNostrum simula- 
tion based on a threshold in the baryon (gas+star) density; it 
was checked that both definitions yield very similar luminos- 
ity functions at various redshifts. Appendix | A4| shows that 
the redshift evolution of the corresponding colours seems 
consistent. 

We choose not to include reddening by dust in the pre- 
dicted colours for the galaxies for two reasons. The first one 
is that our knowledge of dust properties and spatial distri- 
bution with respect to the gas is still somewhat uncertain. 
The accurate modelling of dust attenuation is a complex is- 
sue and we do not want to enter into debates regarding this 
modelling. The second reason is more technical: ray-tracing 
photons involves potentially multiple scattering, and such an 
approach is technically difficult to implement in a very large 
simulation like MareNostrum (see Devriendt et al. ( 2009 1 



and Appendix |A6| for an alternative approach). In the fol- 
lowing, one must keep in mind that the effect of dust on 
galaxies' light is not accounted for in our work. This restric- 
tion should not impact our qualitative findings in terms of 
the influence of filaments, provided dust follows light, which 
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seems a good first order approximation for this simulation 



(Rimes (2009) and Appendix A6l 



We do include, however, IGM absorption on the line of 
sight, which becomes significant in optical bands at z > 2. 



We follow the prescriptions of Madau ( f995) on the hypoth- 
esis of Lya, Ly/3 , Ly7 and LyS line blanketing, induced by 
Hi clouds Poisson-distributed along the line of sight. 

Figure [2] shows examples of galaxies at redshift z = 1.6 
in the MareNostrum catalogue in the I, K and IRAC-8/im 
bands (see Appendix [Bj generated by the spectral synthesis 
described in this section. The consistency of these colours 
with classical models is checked in Appendix |A4| All mag- 
nitudes are in the AB system throughout the paper. Ap- 
pendix [X5] carries a couple of checks on the simulation. 



1 

(5 




2.3 Filaments and the skeleton 

Filaments correspond to the natural framework to charac- 
terise the environments of galaxies on large-scale: within the 
cosmic network, large void regions are surrounded by a fila- 
mentary web linking haloes together. Formally, the skeleton 



(Novikov et al. (20061; Sousbie et al. (2008 1 ) gives a mathe- 
matical definition of the filaments as the locus where, start- 
ing from the filament type saddle points (i.e. those where 
only one eigenvalue of the Hessian is positive), one reaches 
a local maximum of the field by following the gradient. This 
involves solving the equation: 



~dt 



Vp, 



(2) 



for x, where p(x) is the dark matter density field, Vp its gra- 
dient, and x the position. Finding the large-scale structure 
network of filaments involves solving Equation a proce- 
dure which has recently been applied to both data ( |Sousbie 
et al. 2008 \ and simulations ( Caucci et al.|[2008[ ). A sim- 
ilar approach is to classify the different structures (halos, 
filaments, sheets, voids) according to the eigenvalues of the 
Hessian of the density field ( Aragon-Calvo et al.|2007|) or the 
potential (Pogosyan et al. 1998 Hahn et al. 2007 Forero- 



Romero et al. 2009) 



Recently, Sousbie et al. ( 2008 1 introduced a probabilist 



formulation of the skeleton, which amounts to finding the so- 
lution to Equation (J2J at the intersection of the void-patches 
of the density field. The filaments can then be seen as the 
frontiers between the voids. This algorithm has the nice fea- 
ture of constructing a fully connected network of critical 
lines, a crucial feature for this project. The implementation 
of this algorithm on the density field of MareNostrum, with 
a smoothing over 2 /i _1 Mpc yields a hierarchical set of seg- 
ments, where each skeleton segment tracks its connection 
to its neighbouring critical points, together with informa- 
tion relative to the underlying field (density, temperature, 
etc). The result is illustrated on Figure [T] where one can see 
that the skeleton smoothed on these scales traces well the 
large-scale overdense filaments, visible by eye. 



distance [h _1 Mpc] 

Figure 3. Isocontours of the number counts of galaxies in ob- 
served colour versus distance to filament space at redshift z = 1.6. 
A population of redder galaxies is present in the close vicinity of 
the filaments. 
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Figure 4. Average observed colour as a function of distance to 
the closest filament at 2 = 1.6. The error bars, for all the figures, 
represent the one sigma error 011 the average. Galaxies near fila- 
ments tend to be redder on average than galaxies located a few 
megaparsecs away from filaments. 



3 THE INFLUENCE OF FILAMENTS 

The combination of the MareNostrum simulation with spec- 
tral synthesis and the skeleton algorithm allows us to inves- 
tigate the geometric dependence of the spectroscopic prop- 
erties of galaxies on the filamentary environment. 



1 see appendix | A2| for a discussion on the effect of smoothing on 
the results presented here 
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3.1 Gradient of physical properties 

To investigate the influence of filaments on the properties of 
the galaxies, we choose to study the colours in the observer 
frame as a function of the distance to filaments. The locus 
of the filaments (shown in Figure [l]) is computed with the 
above-described skeleton algorithm, using a 256 3 grid and 
smoothing on a — 12 pixels (i.e. 2 ft -1 Mpc). Figure [3] rep- 
resents the distribution of the observed colour G-K (which 
brackets the 4000 A break in the SED at this redshift) as a 
function of distance to large-scale filaments. First, indepen- 
dently to the distance to filaments, the distribution in colour 
shows a bimodality, that will be investigated in further de- 
tails in Section [4]: a distinct population of very red galaxies 
is present. Then it also shows that galaxies tend to be redder 
near filaments. The trend is clearly seen when the distribu- 
tion is averaged (Figure [4|: the G-K colour drops from 2.1 
near filaments to 1.7 at a distance of 5 Mpc. Galax- 
ies exhibit a clear gradient of colour versus the distance to 
filaments. 



3.2 On the influence of nodes and filaments 

The interpretation of this gradient is not straightforward 
since the distance to the skeleton does not only reflect the 
influence of filaments. Indeed, clusters, located at the nodes 
of the skeleton, have been known (e.g. |Goto et al.| ( |2004[ ) 
and references therein) to have a strong influence on the 
properties of the galaxies. Galaxies near filaments are also 
geometrically systematically closer to nodes, and this bias 
could explain the observed gradient. 

The same procedure can be applied to the distance to 
the nodes alone (Figure [5]l . The influence of nodes on the 
colours turns out to be even greater than the effect of the 
distance to filaments; it may thus explain a major part of 
the dependence of the colour with the distance to filaments. 

One can however decrease relatively the contribution of 
the nodes by using a volume average rather than a number 
average: 

i I i 

where pi is the density of galaxies at the position of the i th 
galaxy. Nodes contains most of the galaxies and are there- 
fore over-represented with number averaged weighting, while 
volume weighted means will shift the focus on the filaments, 
which span on much greater scales than clusters. The volume 
averaged colour is plotted against the distance to the fila- 
ments in Figure [5] The colour gradient is strongly damped, 
showing that most of it can be explained by the bias corre- 
sponding to the distance to nodes. 

Even if the influence of nodes is greatly reduced by vol- 
ume averaging, it does not totally vanish; it makes it difficult 
to rule out a weak influence of filaments relative to a resid- 
ual influence of nodes. One would want to know how the 
properties of a galaxy would be modified if its distance to 
filaments was changed while all other parameters, including 
the distance to nodes, are kept unchanged. A way to evaluate 
this is to look at galaxies having the same distance to nodes 
but different distances to filaments. In order to find such 
pairs of galaxies, we proceed in steps (Figure pi. For each 




distance [hr 1 Mpc] 

Figure 5. Circles: average colour versus distance to the closest 
filament at z = 1.6. Triangles: volume averaging. Crosses: dis- 
tance to nodes. The gradient of colour as a function of distance 
to filaments can be mostly explained by the sole influence of the 
nodes and is strongly damped when focusing on the filaments. 




Figure 6. Illustration of the algorithm used to find the filamen- 
tary counterparts. For each galaxy G, we find the closest skeleton 
segment S. We follow the skeleton until the point S" where the 
distance to the node N is sufficiently close to the distance be- 
tween G and N. Then the closest galaxy G' is found and labelled 
as the filament counterpart of G. A galaxy and its filament coun- 
terpart have approximatively the same distance to the node, but 
the filament counterpart is generally much closer to the filament. 
Thus it opens the possibility of unbiased comparison. 



galaxy, we first look at the closest skeleton segment (Gt>S). 
This segment being closer to nodes, we follow the filament 
until we reach a segment with a distance to the node suffi- 
ciently close to the distance between the node and the initial 
galaxy (St>Nt>S'). The closest galaxy to this segment is con- 
sidered to be the filament counterpart of the initial galaxy 
(S' t>G'): it has roughly the same distance to the node, but 
is generally much closer to the filament. The comparison of 
a galaxy and its filament counterpart is therefore of much 
interest to study the influence of the sole filaments. Note 
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distance [h Mpc] 

Figure 7. Difference of observed colour between the galaxies and 
their filament counterparts versus distance to the closest filament. 
Even though a residual 2 sigma shift from zero remains, no overall 
gradient is present. Galaxies far from filaments do not tend to be 
different from their filament counterpart. 



that this construction is not always possible. Indeed if two 
nodes are quite close, and if the intial galaxy is far from 
the filament, the filament linking them is too short to find a 
segment equating the distance to the node. Therefore some 
galaxies do not have filament counterparts and are rejected. 
If we accept only the galaxies which counterpart has a dis- 
tance nodes greater than 90% of the distance of the initial 
galaxy, around two thirds of the galaxies have a filament 
counterpart. Note however that the galaxies rejected are, as 
explained, far from small filaments, while the galaxies close 
to important filaments, which are of interest in our investi- 
gation, should not be affected. Figure [7] shows the difference 
of colour between a galaxy and its counterpart. When this 
procedure is implemented no significant statistical difference 
is found, showing that the influence of the filaments is too 
weak to be detected by this method and that the small gra- 
dient exhibited by the volume averaging procedure can be 
for the most part explained as a residual influence of the 
nodes. 



3.3 Other physical tracers 

The physical properties studied in the previous section are 
the observed colours, since they are easy tracers to observe, 
and are known to reflect well the other properties of the 
galaxies. For the MareNostrum simulation, the same proce- 
dure can however be applied to other physical intrinsic fea- 
tures of galaxies. Figure [8] shows the trend for the rest-frame 
colour UV-I, the specific star-formation rate SFR/M, and 
the galaxy mean stellar age. As for the observed colour, no 
statistically significant gradient remains once the influence 
of nodes is properly removed. 

Figure [9] shows the trend for the metallicity. A gradient 
is still present after volume averaging and can be seen in 
pair comparison. The filaments seem to have an influence 
on the metallicity of the galaxies: the closer to the filaments 
the galaxies are, the more metallic they are. This gradient 
is however fairly small (0.1 dex is the order of magnitude of 
the global change in the metallicity of galaxies within 1 Gyr 
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Figure 10. Distribution of observed colour for several redshifts 
(2 = 1.6,1.75,2.0 and 2.5). The curves are shifted by 2000 for 
clarity. Note the weak bimodality occurring below redshift 2. 



This gradient may nevertheless reflect the enrichment of the 
Warm -Hot Intergalactic Medium (WHIM |C"en fc Ostriker| 
(1999)). Indeed, it is a consequence of the increased galaxy 
density, together with the corresponding (in-)efficiency of 
winds as induced by the subgrid physics. If this result were 
confirmed not to depend critically on the supernovae rate 
recipe, it offers the prospect of directly exploring this com- 
ponent of the IGM via the metallicity of galaxies. It would 
then be of interest to cross correlate this metallicity with 
that of the Orv absorbers (Yoshikawa et al.|2003 l. 



4 BIMODALITY WITHIN CLUSTERS 

In addition to its evolution with the distance to filaments, 
the presence of a bimodality corresponds to an interesting 
feature of the distribution of colours: as seen on Figure [3] a 
population of very red galaxies (G—K > 3) is present. More- 
over, the evolution of the colour distribution with redshift 
shows that this bimodality is appearing around z — 2 (Fig- 
ure 10 1. Let us therefore trace back the positions of galaxies 



responsible for this bimodality in order to explain its origin. 



4.1 Properties of the reddest galaxies 

Figure [TT] shows the distribution of distance to nodes for the 
reddest galaxies (G—K > 3) and for the overall sample. Red 
galaxies are shown to be more clustered near nodes. 

These red galaxies al so e xhibit a very low mass (M 
1O 8 M ). Indeed, as figure 
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at this epoch, see Savaglio et al. ( 2005 1), which is consistent 



shows, the galaxies are either 
blue (with high or low masses) or red with low masses. Such 
small masses associated to very red colours (which in our 
case mean an absence of star formation since we do not in- 
clude internal dust in the modelling of the SEDs) suggest 
that these may be dwarf spheroidals, which origin is still 
poorly known. Grebel et al. ( 2003 1, in a study of dwarf galax- 



with the fact that the other properties, which are indirectly 
related to metallicity, do not exhibit any significant gradient. 



ies in the Local Group, suggested that the absence of star- 
formation in such objects is due to externally induced gas 
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Figure 8. Same as Figure [5] for UV-I rest-frame colour (left), specific star formation rate (middle) and mean light-weighted stellar age 
(right). As expected, the same trends are found for these physical tracers as for the intrinsic colours. 




distance [h~ 1 Mpc] distance [h 1 Mpc] 



Figure 9. Same as Figure[5] (left) and Figure[7] (right) for the metallicity. Metallicity shows a clear gradient as a fuction of distance to 
filaments. This gradient reflects the inhomogcncity of the WHIM on these scales for the MareNostrum subgrid physics. 



loss: ram-pressure stripping would be responsible for lack of 



ISM gas in these small galaxies. Mayer et al. (2007 1 showed 



that stripping in the clusters can indeed increase the mass- 
to-light ratio of dwarf galaxies. It is possibly the case here 
too, as demonstrated in Sect. |4.2] and Figure [13] In contrast 
to what is observed at low redshift, there is no massive red 
galaxy in the simulation at z — 1.6. [^] 



2 Low mass galaxies should be considered with caution in such a 
simulation. For instance low mass galaxies are expected to be ar- 
tificially blue as they only have a recent accretion history, whereas 
more massive galaxies are spectrally better resolved. However this 
effect would not explain the fact that the galaxies studied are red- 
der. 



4.2 Dynamical scenario 

The dependence of the spectro-photometric properties with 
the distance to filaments studied in the previous sections re- 
flects the process of galaxy formation and its connection 
with the global flow within the large-scale structure. In- 
deed, it has been demonstrated ( Aubert et al.|2004 Sousbie 
|et al.|2008[ ) that filaments are fed by surrounding voids and 
mark the lanes of galactic infall towards the clusters. Young 
galaxies form across the whole filamentary network, but are 
rapidly collected along the more busy subnet of denser fila- 
ments. 



A close-up look at a specific halo (see Figure 13 1 con- 



firms the presence of small, red galaxies (which do not form 
stars) embedded in the halo. These galaxies are already 
accreted at z = 2. At this range of redshift, haloes con- 
tain big (M > lO lo A/0), blue, central galaxies and small 
(M < 1O 9 M0), red galaxies which are being stripped of 
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Figure 13. View of a small region of the simulation centred on a large cluster. The traces represent the position of the galaxies for 
several timestcps between 2 = 1.9 and 2 = 1.6. Each galaxy is represented by a point whose size depends on its mass and whose colour 
represents either the SFR/M* (left) or the UV-I rest-frame colour (right). The top figures represent a 4 h _1 Mpc cube and the bottom 
ones are a zoom on the bottom right part and span over 1 h _1 Mpc. The grey background encodes the temperature of the gas, showing 
the extent of the hot gas bubble. Note the reddening of the small galaxies entering the cluster (bottom right), in parallel to the reduced 
star formation (bottom ieft). 



their gas and swallowed by the central galaxy. The fate of 
these small galaxies can be investigated by looking to an- 
other population: the galaxies that are currently accreting 
into the halos at this epoch, although they have a larger 
mass (1O 9 M < M < 10 10 M o ). These intermediate galaxies 
shows a peculiar behaviour that could explain the observed 
increase of the average colour near nodes: when approach- 
ing the halo, they stop forming stars and therefore begin to 
passively redden. This quench could be triggered by the lack 
of cold stream reaching the inner core o f dark halos (|Ocvirk| 
et al.|2008| o r by stripping of their gas ( |Moore et al.|2000[ ). 
Ocvirk et al. ( 2008 1 show that more massive haloes prevent 



cold stream from feeding the inner halo; here it is found that 
these halos display redder satellites. The accreted galaxies 



tend therefore to be red and have a low SFR, leading to 
the previously observed bimodality. They would then merge 
with the central galaxies, and contribute to the increase in 
mass (dry merging) of the massive old galaxies observed at 
low redshift. 

In this scenario, the spectroscopic properties of a galaxy 
should be at least in part determined by the physical con- 
ditions of the intra-cluster medium. Figures [14] and [15] are 
consistent with this scenario. Indeed the average colour in- 
creases with the pressure, as expected in the case where the 
reddening is induced by star formation quenching via ram- 
pressure stripping. The detailed 2D distribution shows the 
presence of a population of very red galaxies for high pres- 
sures. This galaxies could correspond to the galaxies which 



On the filamentary environment of galaxies 9 



0.25 



0.20 



£0.15 



0.10 



0.05 



red galaxies (G-K>3) 

H III 



i ' i ■ i 1 1 1 1 1 1 1 1 1 1 




all galaxies 



,1,1,1,1,1,1 



2 4 6 8 

distance to nodes [h~ 1 Mpc] 



Figure 11. PDF of the distance to nodes for all galaxies (black) 
and reddest (G — K > 3) galaxies (red). The reddest galaxies are 
clustered near the nodes. 
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Figure 15. Distribution of galaxies in observed colour (G-K) 
pressure space. The increase of the G-K colour with pressure is 
due to a distinct population. 
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Figure 12. Isocontours of the density of galaxies in mass versus 
colour space, for z = 1.6. Red galaxies are typically less massive. 
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Figure 14. Mean observed colour, G-K, as a function of pressure 
of the galaxy. Galaxies are redder in high pressure regions. 
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Figure 16. Pressure of the galaxies versus distance to nodes. 
High pressures (likely to strip low mass galaxies of their gas) are 
only found near massive clusters. 



have been stripped of their gas when their entered the clus- 
ter, and will then become red, small galaxies, before they 
eventually merge with the core galaxies. 

The typical pressure of transition is around 10~ 14 Pa 
and corresponds to conditions found only near nodes (see 
Figure|16[). Theoret ical considerations and simulations ( |Fu-| 
jita & Nagashima 1999 Roediger 2009 1 show that ram- 



pressure stripping is expected to be efficient for pressure 
over 10~ 13 Pa. This apparent discrepancy may be explained 
by noticing that we are considering different objects. Their 
value corresponds the pressure required to strip a typical 
spiral galaxy at low redshift and could thus be lower for the 
smaller high redshift galaxies we are considering. 
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4.3 Tentative link with observations 

It is clearly beyond the scope of this paper to carry a di- 
rect comparison with the currently available data. Indeed 
the main difficulty in comparing the results found in this 
simulation with existing observations lies in the range of 
redshifts and the geometry of the surveys. Building a fully 
connected skeleton is non trivial for pencil shape volumes, as 



edge effects become important (as discussed in Sousbie ct al 
( 2008 1). Moreover, in this paper, the skeleton was computed 
from the dark matter distribution and one would need to 
calibrate the bias involved in using light instead of mass. 

Finally, the simulation was stopped atzw 1.5, and the 
bimodality seems to appear around z as 2 (see Figure 10 1. 



Observations at such high redshift are uncommon and diffi- 
cult (|van Dokkum et al.|2006[ |Daddi et al.|2005 l. However, 



as mentioned in Section [1] similar trends exist for the same 
range of redshift: for example, red galaxies are observed to 
be more clustered ( [Daddi et aL]|2003[ ). 

These gradients can also be linked to slightly lower red- 
shift observations. For redshift z f» 1, there exists several 
surveys that have studied the role of the environment, e.g. 
GOO DS ( |Elbaz et al.p007||, VVDS ( |Scodeggio fc Vergani 
2009[ ) or DEEP2 ( |Cooper et al.|2006[ ). The existence of a bi- 
modality, with a blue and a red sequence, is known to appear 
at this redshift ( Nuijten et al.|2005 1 and can be observed on 
several properties of galaxies at low redshift ( Mouhcine et al. 
[20071). 



To better understand the link with observations, one has to 
understand the robustness of the results with respect to ob- 
servational uncertainties and other biases. The nature of this 
virtual data set allows us to compute directly the proper- 
ties (skeleton, distance to filaments, . . . ) in 3D, without any 
problems of distance determination. For observational ap- 
plications, a natural question arises: are these results robust 
with respect to distance uncertainty, that can come e.g. from 
the use of photometric redshifts ( jCooper et al.|2005[ )? In or- 
der to asses this robustness, the simulation cube is projected 
along one of its axis. The 50 h~ l Mpc can be thought of as 
representing the uncertainty on distance (roughly Sz — 0.005 
at z — 1.6). A 2D skeleton is computed on the resulting 
density field, with the same N-dimensional algorithm. The 
observed colour is then compared to the 2D distance to the 



closest filament (see Figure 17 1. The main features found in 
the 3D investigation are still weakly present after projection: 
a distinct population can be seen near nodes and leads to 
a enhanced averaged colour, but filaments do not seem to 
have an effect on the properties of the galaxies. This effect 
is however far more subtle than in 3D, given the projection 
effect. 

Finally, as an alternative to the skeleton as a tracer of 
filaments let us briefly implement a few more commonly used 
tracers (as mentioned in the introduction) on the MareNos- 
trum simulation: the 2D distance to the 5 th neighbour 
(|Cooper et al.|2005 \ and the galaxy number density. Figure 



18 shows that a bimodality is still present with the 5 neigh- 
bour probe, but it is much less contrasted and much more 
localised. The 2D distance to the 5 th neighbour is therefore a 
less sensitive probe of the anisotropic cosmic environment. 
In contrast, provided spectroscopic redshifts are available 
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Figure 17. Same as Figure[5]for the 2D skeleton. A weak trend 
is still present but much less significant than in 3D. 
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(say with the LSST, Claver et al. 1 20041), the 3D skeleton 
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Figure 18. Isocontours of the number count of galaxies in ob- 
served colour versus distance to the 5 neighbour for the pro- 
jected data at z = 1.6. Note the larger number of contours used 
to catch the low contrast bimodality, compared to Figure [3] Note 
also the difference of scale on the x axis: the 5 th neighbour does 
not probe large-scale structures when no cut on galactic mass is 
applied. 



should allow us both to probe the large scale structures and 
mark the neighbourhood of clusters in detail. 

Finally, figure [19] shows that the bimodality is not 
present anymore when the galaxy number density is used, 
in 2 or 3D. 



5 CONCLUSION & DISCUSSIONS 

The Cosmic Web is a key feature of the organisation of 
galaxies on large scales. This paper investigated the influ- 
ence of this filamentary environment on the spectroscopic 
properties of galaxies, using the MareNostrum simulation 
which was postprocessed using stellar population synthesis. 
The cosmic web was traced with the skeleton algorithm, and 
has proven here to be a very effective mean of probing the 
anisotropy of the large-scale structures. 

We found gradients of spectroscopic properties of galax- 
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Figure 19. Isocontours of the number count of galaxies in ob- 
served colour versus number of galaxies within a sphere of 1.5 
h~ x Mpc for the projected data (top) and the 3D data (bottom). 
No significant trend is seen with these classical tracers. 



ies with the distance to filaments, but demonstrated that 
they can be explained by the fact that the distance to fila- 
ments is biased by the distance to the nodes of the network 
(group or clusters of galaxies). Two procedures were intro- 
duced to remove the influence of these nodes: (i) volume- 
averaging this distance decreases the influence of the com- 
pact, dense regions and focus on wider structures, such as 
filaments; (ii) pair comparison which seeks the filamentary 
counterpart of each galaxy. Both methods show that the 
influence of the filaments alone is negligible compared to in- 
fluence of clusters. A bimodality in colour was also found to 
occur below redshift z ~ 2 and its origin was investigated. It 
is due to a population of red, small galaxies (~ 1O 8 M0) ac- 
creting on the nodes of the Cosmic Web, while more massive 
objects, M > 1O 8 M0, are mostly unaffected. These galaxies 
have their star formation quenched while they approach the 
clusters. It remains to be confirmed that this stripping pro- 
cess is not amplified by a lack-of-resolution effect, since (i) it 
involves amongst the smallest (virtual) galaxies in the simu- 
lation, and (ii) it seem to create a tension with observations 
at redshift zero ( |Kimm et aL|2009| ). 

These findings suggest that the large-scale filaments are 
only dynamical features of the density field, reflecting the 
flow of galaxies accreting on clusters; the conditions in the 
filaments are not dramatic enough to influence strongly the 
properties of the galaxies it encompases, unlike the intra 
(proto) cluster medium, which seems able to strip down the 
ISM of the incoming low mass galaxies. Appendix |A3| shows 
that this statement remains valid when the study is limited 
to the most important filaments. 

Finally, we did find a weak metallicity gradient away 



from the filaments which could reflect the large-scale in- 
homogeneity in the distribution of metals (the so called 
WHIM). 

One could have imagined that even if the large-scale 
filamentary network were purely a tracer of the large-scale 
dynamics, galaxies within the large-scale filaments should 
be redder and older, since they would have joined the cos- 
mic super highway earlier on average when compared to field 
galaxies. This effect is not seen at those redshifts, as (i) older 
galaxies continue to accrete new cold gas on small scales and 
form stars, hence remain blue, (ii) some field galaxies contin- 
uously join the large filamentary network, and (iii) on large- 
scales, a significant fraction of the matter is accreted more 
or less radially onto the nodes, while filaments are collecting 



left-overs from this accretion more indirectly (Sousbie et al 



2008) 



Recently (Keres et al. 2005 Ocvirk et al. 2008 Dekel 



et al. 20091, it was emphasized in steps that anisotropic 



metal-rich cold stream accretion regulates the inflow of cold 
gas towards the inner regions of the most massive galaxies of 
the high redshift universe (z > 2). It was then conjectured 
that this process could explain the observed bimodality of 
spectroscopic galactic properties at lower redshift. In this 
paper, we have shown that the geometric distribution of the 
colour of galaxies is not sensitive to the detailed large-scale 
filamentary network, but only to its nodes. This apparent 
paradox may be lifted when noting that the self-regulating 
anisotropic filamentary accretion occurs on much smaller 
scales, and was quantified for the central galaxies of the 
simulation which are sitting at the nodes of the network; 
in contrast, when considering the full galactic population, 
a typically low mass galaxy is not transformed by its en- 
counter with the different physical condition of the weakly 
overdense intergalactic medium within the cosmic web, un- 
less it falls into the intra cluster medium of a large node. 
In other words, massive galaxies feel the small scale filaments 
(cold streams) feeding them at nodes; low mass galaxies are 
not spectroscopically changed while entering the large-scale 
filaments. The mesoscopic (below a Mpc scale) filamentary 
feature of the cosmic network may geometrically solve the 
self-regulating process of galactic accretion, but we have 
demonstrated here that its large-scale counterpart does not 
seem to directly affect the colours of galaxies. 

In Sousbie et al. (20081, the effect of redshift distortion 



is partially addressed, and the corresponding algorithm is 
now being extended to discrete surveys via a Delaunay tes- 
sellation (which are therefore not sensitive to edge effects). 
This, as argued in Section [4. 3| is a critical step towards per- 
forming a similar 3D analysis on real data, which as we have 
shown is essential to quantify these gradients, as in projec- 



tion, the information is lost (see Figure 17 1. In particular, it 
would be of great interest to carry out these measurements 
on a DEEP-2/VVDS/z-COSMOS-like survey as well as to 
bring the simulation down to lower redshift and reach a time 
in cosmic history when upcoming large observational sur- 



veys (e.g. LSST Claver et al. (20041, BOSS Schlegel (20071) 



overlap statistically with the predictions of the simulation. 
It would also be worth investigating how sensitive some of 
our findings are w.r.t. the detailed chosen subgrid physics 
by probing alternative recipes and running higher spatial 
resolution simulations. 

The catalogues produced for this investigation (spectro- 
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scopical properties of the MareNostrum galaxies and its dark 
matter skeletons) are available online as discussed briefly in 
Appendix [B] 
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APPENDIX A: SELF-CONSISTENCY CHECKS 

In this section, we present a few tests that were performed 
to check the presence of bias and the influence of the pa- 
rameters of the skeleton and the MareNostrum simulation. 

Al Decorrelating galaxies with their environment 

In order to check if our method induces artificial correlations 
of the galaxies properties with their environment, we shuffle 
the properties of the galaxies: for each galaxy, we change 
its colour to the colour of another randomly chosen galaxy, 
keeping its position unaffected. The results obtained on the 
shuffled sample, corresponding to Figure [5j are presented in 
Figure |A1| As expected, the evolution of colour with the 
distance to filaments or nodes is totally erased. 
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Figure A2. Left: Galaxies (red) and skeleton (blue) for a smoothing length of 8 (top) and 16 (bottom) pixels. Increasing the smoothing 
length selects the main features of the cosmic web. Middle: Same as Figure[5]for a smoothing of 8 (top) and 16 (bottom) pixels. Right: 
Same as Figure[7]for a smoothing of 8 (top) and 16 (bottom) pixels. 
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Figure Al. Same as Figure[5] (top) and Figure [7] (bottom) when 
galaxies properties are shuffled. As expected, the evolution is can- 
celed. 



A2 The influence of smoothing 

To compute the skeleton, the density field needs to be 
smoothed to insure sufficient differentiability. The smooth- 
ing length is a free parameter, allowing one to probe differ- 
ent scales. All the previous results have been obtained with 
a Gaussian smoothing over a =12 pixels ( with a 256 3 grid), 
which corresponds to 2 h~ 



A2 



shows the re- 



Mpc. Figure 

suits for a smoothing length of 8 and 16 pixels. Increasing 
the smoothing allows to select only the biggest features of 
the cosmic web. Thus, it changes the influence of the nodes 
in two different ways. First it increases the average colour 
near nodes, the reddest galaxies being located in the biggest 
clusters. Then it increases the influence scale of the nodes: 
galaxies properties can be influenced by the main clusters 
even at distances of several megaparsecs. 

Nevertheless, the influence of the filaments is still van- 
ishing when the influence of the clusters is correctly re- 
moved. Even the biggest filaments do not have a direct effect 
on the properties of galaxies. 



A3 Spurious filaments 

The fact that that some filaments of the skeleton could be 
spurious and would not correspond to any physical filaments 
could be responsible for the lack of dependence of spectro- 
scopic properties with the distance to skeleton. This would 
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Figure A4. Evolution of the distribution of the rest-frame colour. 
Black lines: maximum and secondary maximum of the distribu- 
tion. Green line: PEGASE model. 



lead to a dilution of the dependence, since the galaxies near 
spurious filaments would not be correctly taken into account. 
In order to check that the skeleton algorithm is not intro- 
ducing such spurious filaments, we remove the less physical 
filaments. Considering whole filaments, i. e. a set of contigu- 
ous segments between two given nodes, ensure us that this 
procedure will not depend on the distance to nodes and will 
not introduce any bias. The selection criterion is based on 
the density of the underlying field, averaged over the fila- 
ment. We choose a threshold such that up to 50% of the 
skeleton is removed. The result of the pair comparison pro- 
cedure is given in Figure |A3| It shows that even when the 
less significant filaments produced by the skeleton algorithm 
are removed, the colour of the galaxies does not seem to de- 
pend on the distance to filaments. The result presented in 
the main text is therefore robust and cannot be explained 
as an artifact of the skeleton algorithm. 



A4 A toy model for the observed colours 

Deriving colours for objects in the simulation is a long and 
complex process, initiated from a model for the primor- 
dial density fluctuations, and involving huge computing re- 



sources to grow the structures. As mentioned in section [2~T) 
various approximations and recipes are required to obtain 
the end-products of interest here, namely (virtual) galaxies. 

To check that the colours presented here make sense 
with respect to what we know of galaxy colours today, we 
can either compare them to observations or to previous mod- 
els. For simplicity, and to avoid potential biases inherent to 
observations done at various redshifts, we choose to test our 
results against a reasonably simple model of galaxy forma- 
tion. Another reason to do so is because we did not include 
dust in the SED modelling of the simulation. A direct com- 
parison with observations is therefore hazardous. 

As a basis for this comparison, we use an idealised sce- 
nario with a smooth star-formation history leading to the 
average colours of local late spiral Sd galaxies. Such a sce- 



nario is presented in Fioc & Rocca-Volmerange ( 1999 1 or 
|Le Borgne et a£\ ( |2004[ ). It involves the infall of primordial 
gas from a reservoir onto a potential well at a rate propor- 
tional to exp(— £/t)/t with r = 6 Gyr. As the accreted gas 
cools, stars begin to form with a Schmidt Law, at a rate pro- 
portional to the gas density SFR= (14Gyr) -1 p ga s. Modelled 
with the code PEGASE. 2, the metallicity of such a galaxy 
evolves consistently with the yields from supernovae, and 
many physical properties are monitored, as well as predicted 
spectra and colours. These two time-scales are the only pa- 
rameters that were tuned to match the average colours of 
local star-forming galaxies with the hypothesis that the red- 
shift of formation of the first stars is z = 10. In practice, 
choosing z > 4 produces an almost equally satisfying match 
for the observed colours. 

This scenario, combined with others, is also successful 
in reproducing the galaxy counts in optical and NIR bands 
( Fioc fc Rocca-Volmerange||1999 1. Although such successes 
are appealing, they do not prove that the scenario is cor- 
rect and it must be taken with caution. Still, we venture 
into the comparison of the MareNostrum simulated colours 
with the evolving colours derived from such an idealised sce- 
nario. We compare in Figure |A4| the evolving distribution 
of the UV-I rest-frame colour (which encompasses the 4000 
A break) for galaxies in the simulation with this so-called 
"monolithic" scenario. For consistency with the simulation, 
we do not include dust in the model for the Sd spiral. We 
find a remarkably good agreement between the colour peak 
of the blue sequence and the colour of the scenario at every 
redshift between z — 1.6 and z — 3, in a range where the 
refinement of the grid used for the simulation is comparable 
from one bound to the other. It suggests that the bulk of 
the galaxies might follow the path of the idealised scenario 
and end-up in local late-type spirals. Such an agreement is 
comforting and suggests that the colours derived from the 
simulation are very reasonable in this redshift range. We 
must stress, however, that this comparison does not vali- 
date the monolithic scenario which was mainly tuned to re- 
produce z — colours: the match of its colours with distant 
galaxies is much more difficult to check because of the dif- 
ficulty to identify the progenitors of local spirals at higher 
redshift. And of course, this comparison cannot be used to 
argue that this cosmological simulation can be reduced to a 
monolithic collapse, which is contrary in nature to the hi- 
erarchical paradigm. Nevertheless, we might learn from this 
exercise that on average, the evolution of the population of 
blue galaxies in MareNostrum is driven by an average accre- 
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Figure A5. Comparison between the star formation rate in the 
MareNostrum simulation and in observed surveys popk ins Hi 
|Beacom"||2006| l . The virtual SFRD seems globally consistent, if 
somewhat high. 



u 
a. 



cn 
D 

E 



O 








i i i i i i 


-1 




O Calzetti, E(B-V)=0.22 j 
A Bouwens et ol (2007) 


-2 






-3 






-4 






-5 


; 


z = 4.0 | i 

: 


-6 




. i ... i ... i ... i ... i . 




-16 


-18 -20 -22 -24 -26 



1400W 

Figure A6. the UV (1400 A) rest frame galaxy luminosity func- 
tion (LF) measured in the Mare Nostrum simulation at redshift 
4. Red diamond LFs stand for galaxies uniformly extinguished 
with a Calzetti (2001) law. Blue triangles correspond to data 
gathered by Bouwens and collaborators | |Bouwens et al.j|2007| ). 
Vertical dashed lines indicate mass resolution, horizontal dashed 
lines volume resolution limits. 



tion rate following the natural law used above, and that on 
average the Schmidt law described above is representative 
of the global star-formation activity. Still, the exact values 
of the time-scales used for our scenario should not be taken 
for granted: a significant degeneracy exists between them, 
not to mention the hypothesis made on the universality of 
the IMF. 

The main outcome of this test is that the rest-frame 
colours that we produce seem coherent, and are consistent, 
if we extrapolate them to z = 0, with the colours of local 
star-forming galaxies. 

A5 MareNostrum simulation sanity check 

Let us carry a couple of checks on the features of the 
MareNostrum simulation and its postprocessing using clas- 
sical probes. Figure |A"5] compares the star formation rate in 
the MareNostrum simulation and in the [Hopkins fc Beacom| 
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Figure A7. The effect of dust on G-K colours as a function 
of stellar mass at redshift 4. Dust is modeled from the metallic- 
ity of a fair sub-sample of the underlying galaxy following the 
prescription of |Devriendt et ah] ( |2009fr . Note that the low mass 
objects responsible for the bimodality described in this paper are 
not reddened. 



( 2006 1 compilation. The virtual cosmic SFRD seems roughly 
consistent with the observed one (which is corrected for dust 
reddening), although on the high end. Figure A6 



compares 



the UV (1400 A) rest frame galaxy luminosity function (LF) 
for galaxies measured in the Mare Nostrum simulation at 



redshift 4 uniformly extinguished with a Calzetti (2001 1 law 
of E(B — V) = 0.22, to that observed by Bouwens et al. 
(20071. Our simulated LFs at z~4 are in good agreement 



with the available UV data (see Devriendt et al. (2009) for 
a more detailed comparison). The apparent discrepancy be- 
tween the perfect agreement of the UV LFs and the slight 
overprediction of the SFRD can be explained by several fac- 
tors: the z = 4 data points for the SFRD from |Ouchi et al~] 



pooi i are corrected with a dust extinction measured on in- 
dividual objects, spanning a range of E(B-V) ~ 0.0 to 0.4. 
Moreover, other uncertain factors are involved in the com- 
parison such as the correction for the IGM absorption or 
the model-dependent law to convert from UV luminosity to 
SFR. 



A6 Tentative modeling of the effect of dust 

Let us finally investigate the effect of dust on the G-K 
colours computed in this paper. Following |Devriendt et ah] 
(2009), we compute the internal absorption of each galaxy 
using the metallicity of the gas within that galaxy as a proxy. 
Figure [A7] shows the difference in G-K colours between the 
extinct and the non-extinct galaxies of the MareNostrum 
simulation at z — 4 as a function of galactic stellar mass. 
It reaches ~ 1 magnitude at ~ lO^M©, but remains quite 
small for the low mass satellites responsible for the observed 
bimodality. 

Figure |A8| displays the same colour gradient as Figure [5] 
(and Figure [7| at redshift 4 for reddened G-K colours; the 
corresponding unobscured colour gradients are shown as 
dashed lines. Taking into account the dust (i) produces a 
global shift of the colours, and (ii) increases somewhat the 
colour gradient as a function of the distance to nodes, as 



massive (more obscured, see Figure A7l galaxies are statis- 
tically more present at the nodes. Nevertheless the effect of 
dust on filaments remains much lower than the influence of 
nodes. 
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Name 



Reference 



Mean wavelength 



FUV Galex 

G 

R 

I CFHT 12K 
SDSS-z 
Johnson K 

SPITZER IRAC channel 1 IRAC-3.6/im 
SPITZER IRAC channel 4 IRAC-8/im 



Steidcl & Hamilton ( 1993 ) 
Stcidcl & Hamilton ( 1993 ) 
Le Fevre et al.|l|20o4 



Morrissey| \200b\ 



Fukuf 


;ita ct al. ( 1996 


Johnson et al. 


(1966) 


Fazio 


{2004} 


Fazio 





1520 A 
4810 A 
6980 A 
8130 A 
8960 A 
21950 A 
35610 A 
79580 A 



Table Bl. Avalaible filters. 



Name 


Description 


Unit 


First 


Last 


xpos 


X position 


Mpc/h 


0.201926 


1.3122 


ypos 


Y position 


Mpc/h 


16.6967 


45.9471 


zpos 


Z position 


Mpc/h 


47.9065 


43.5113 


rcdshift 


Redshift 


- 


1.56506 


1.56506 


zJ 


Luminosity-averaged stellar metallicity 


dex 


-1.69139 


-2.97829 


z_m 


TV T J _i 11 , 1 1 ■ ■ i 

Mass-averaged stellar metallicity 


dex 


-1.78234 


-3.03794 


age-1 


Luminosity-averaged stellar age 


Myr 


931.873 


1429.32 


age_m 


Mass-averaged stellar age 


Myr 


1763.61 


1640.19 


FUV 


Observed AB magnitude in the FUV band 


mag 


28.0128 


47.3694 


G 


Observed AB magnitude in the G band 


mag 


20.9583 


36.5068 


K 


Observed AB magnitude in the R band 


mag 


20.9726 


35.7524 


1 


Observed AB magnitude in the I band 


mag 


20.8269 


35.2413 


z 


Observed AB magnitude in the z band 


mag 


20.6424 


34.8702 


K 


Observed AB magnitude in the K band 


mag 


19.3598 


33.2805 


IRAC3p6 


Observed AB magnitude in the IRAC 3.6 fim band 


mag 


19.1845 


33.252 


IRAC8 


Observed AB magnitude in the IRAC 8 fim band 


mag 


20.0335 


34.2273 


FUV_restframe 


Absolute (rest-frame) AB magnitude in the FUV band 


mag 


-25.3455 


-9.4258 


G_rest frame 


Absolute AB magnitude in the G band 


mag 


-26.5264 


-12.6057 


R_rcstframe 


Absolute AB magnitude in the R band 


mag 


-23.829 


-12.9374 


I_restframe 


Absolute AB magnitude in the I band 


mag 


-26.9103 


-13.0224 


z_restframe 


Absolute AB magnitude in the z band 


mag 


-27.0163 


-13.0619 


K_restframe 


Absolute AB magnitude in the K band 


mag 


-36.8808 


-12.6582 


IRAC3p6_rest frame 


Absolute AB magnitude in the IRAC 3.6 fim band 


mag 


-26.1501 


-11.8306 


IRAC8_restframe 


Absolute AB magnitude in the IRAC 8 fim band 


mag 


-24.7171 


-10.2381 


distance_skel 


Distance to the skeleton 


Mpc/h 


0.371031 


3.73985 


distance_node 


Distance to the nodes of the skeleton 


Mpc/h 


0.811432 


5.69763 


sfr 


Star Formation Rate 


M Q /yr 


95.8817 





stellar _mass 


Total stellar mass 


M Q 


5.11478e+ll 


9.93592e+06 



Table B2. Content of the catalogues. The first and last galaxies of the catalogue at z = 1.57 are given as reference. 



APPENDIX B: PHOTOMETRIC CATALOGUES 

The data used in this study are made publicly available 
online: http : //www. iap.fr/users/pichon/MareNostrum/ 
catalogues. The catalogues contain the properties of the 
galaxies (position, colour, metallicity, age, SFR, stellar 
mass) and their environment (distance to skeleton and to 
nodes, corresponding densities) for redshifts 1.57, 1.8, 2.1, 
2.51, 3.01, 3.53 and 3.95. The number of galaxies in each 
catalogue is respectively 97563, 103589, 111184, 119978, 
124642, 119612 and 103187. See Table|B2]for a detailed list 
of the properties. The observed and the rest-frame colours 
are given by the magnitudes in the 8 filters described in 
Table [Blj Other filters could be implemented upon request. 
For the sake of simplicity, the skeleton are saved as a set 
of coordinates together with the mean dark matter density. 
The node catalogue follows the same prescription. Here the 
catalogues are distributed as VOtable^ FITS and ASCII. 



http: / / www.ivoa.nest/Documents/latest / VOT.html 



These could be of use to anyone intending to compare 
this large-scale simulation to e.g. observations at high 
redshifts. The skeletons are also available there for further 
investigation. 



APPENDIX C: SUBGRID PHYSICS 

The simulation analyzed in this paper relies on subgrid 
physics which we now briefly summarize. 



CI Star formation rate 

Let us recall the method for implementing star formation 



in RAMSES as described in Rasera fc Teyssier (20061. It is 



based on a phenomenological approach adopted in many, if 
not all, cosmological studies. It relies on observations of the 
scaling law between star formation rate and gas densities 
of local galaxies ( |Kennic utt 1998). Basically, one considers 
that star formation proceeds at a given time scale, written 
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Figure A8. Same as Figure [5] (top) and Figure [7] (bottom) at 
higher redshift (z=4) with (plain line) and without (dashed line) 
dust absorpion. The net effect of dust is a global offset of the 
colours. The colour gradient for the nodes is slightly stronger with 
dust, which reflects the presence of massive dusty galaxies at the 
nodes, but does not modify our conclusions: filaments alone do 
not have a significant influence on the (dusty) colours of galaxies 
(bottom panel). 



where m* is the final mass of the star particle, t?sn is the 
fraction of mass in supernovae ejecta per solar mass of stars 
formed (jjsn = 0.1 for the assumed Salpeter IMF) while 
rjw is the mass loading factor that determines the gas mass 
entrained by the supernovae ejecta (which is set here to 1, 
see Dubois & Teyssier (2008)). We assume that these debris 



are distributed according to a Sedov blast wave solution with 
a maximum speed given by 



MSN 



\Jl + Vw/vs-t 



(C4) 



where «sn is the typical velocity corresponding to the kinetic 
energy released in one single supernova explosion (wsn — 
3200 km. s -1 ). The energy released to the gas by the debris 
is 



Ed = ??SN 



M SN 



(C5) 



where Msn and -Esn are respectively the typical progeni- 
tor mass and energy of an exploding type II supernova (i.e. 
M SN = 10 M Q and E S n = 10 51 erg). 



here t* , in regions where one or several physical criteria are 
fulfilled. We adopt a simple scheme to turn gas mass into star 
particles, adding a source term in the continuity equation 



(Dp\ 
\DtJ* 



-p/t* , if p > pa = 0.1cm 
, otherwise, 



(CI) 



where the star formation time scale i» is proportional to the 
local free-fall time 



U = to 



-1/2 



with to — 2Gyr. 



(C2) 



This choice of to corresponds to the 5 % efficiency mentioned 
in the main text. 



C2 Supernovae feedback 



Following closely Dubois & Teyssier (20081, we define the 

(C3) 



mass vanished in the star formation process as 



(Am 9 ) SF =m«(l + rjsN + Vw) , 



